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Abstract 

The process of charged Higgs production in association with a top quark at the LHC has been 
calculated at the complete NLO electroweak level both in a Two Higgs Doublets Model and in 
the Minimal Supersymmetric Standard Model, assuming a mSUGRA breaking scheme. We have 
numerically explored the size of the one-loop corrections in two typical supersymmetric scenarios, 
with particular attention to the tan/3 dependence, and we have found that they remain pertur- 
batively small but possibly sizable, reaching a 20% limit for extreme values of tan /?, when the 
complete set of Feynman diagrams is taken into account. 

PACS numbers: 
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I. INTRODUCTION 



The processes of production of a charged Higgs boson will be extensively exploited to 
search for new physics beyond the standard model at the LHC. Most extensions of the 
Standard Model (SM), such as two-Higgs-doublet models (2HDM) or the Minimal Super- 
symmetric Standard Model (MSSM), enlarge the minimal SM Higgs sector predicting the 
existence of charged Higgs particle(s). Since the discovery of a charged Higgs boson would 
be a distinctive signature of new physics, an exhaustive comprehension of its production 
mechanism appears to be mandatory. 

Depending on the charged Higgs boson mass, different production mechanisms are dom- 
inant: if mf^+ < rrit — mi, the main source of charged Higgs is the tt production and the 
subsequent decay of the top t H^h hj, while for a heavier charged Higgs boson the dom- 

n n 

inant process is the associated production with heavy quarks [2|-[12|. Also the associated 
production with W gauge boson has been analysed 3]- [3], but this process is suppressed 
with respect to the other two mechanisms of production. 

We will focus our analysis on the associated production with a top quark, which is also 
an important mechanism of top production and should be considered in the analysis of 
single top production at the LHC [l^. At the lowest perturbative order, it is well known 
that this process is particularly sensitive to the value of the parameter tan/5, i.e. the 
ratio of the neutral Higgs vacuum expectation values V2/V1., which appears in the Yukawa 
coupling thH . Being proportional to mbtan/5, the coupling is enhanced for large values of 
tan (3 and this enhancement allows a direct check of the 2HDM structure of the model, not 
necessarily involving SUSY. Supersymmetric corrections, on the other hand, can be only 
investigated looking at the loop structure of the process. 

Due to its relevance, the process of production of charged Higgs in association with a top 
quark has been extensively studied at higher orders and many important results have been 
obtained. The NLO corrections in QCD and SUSY QCD in the five-flavour scheme {i.e. 
including the bottom quark as a parton of the sea) have been computed in Refs. [3]- [ill ^^"^ 
the same corrections in the four-flavour scheme, together with a comparison of the results 



in the two schemes have been computed in Ref. 



12|. As a general feature, while QCD 



corrections are generally found to be large and positive and nearly independent of tan/3, 
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SUSY corrections appear to be sizable and negative for large tan [3. For what concerns the 
NLO electroweak (EW) contribution, the subset of Yukawa SUSY EW corrections has been 
computed in Refs. Qj-Q] in both the five- and four- flavour schemes: all of these papers 
assume that the Yukawa part of the correction is the leading one and they get some large 
one-loop contributions for "extreme" values of tan (3. 

Given the possible relevance of the considered process, which might require a more accu- 
rate prediction, we have performed in this paper a complete NLO MSSM EW calculation. 
This includes all the EW diagrams that were neglected in the previous analyses and also the 
total QED radiation, that has never been computed for this process and whose effects might 
be a priori relevant, as we know from previous recent calculations of our group 17|, ll8|. 

The paper is organized as follows. Sect. [Illwill be devoted to a description of the shape 
and of the basic properties of the parton level amplitudes for hg tH~ at Born and at one- 
loop level. A rigorous treatment of QED radiation has been performed to obtain reliable 
values for the observables and will be described in Sect. Ill CI In Sect. IIIII the numerical 
one-loop effects on the production rates and distributions for a couple of meaningful SUSY 
benchmark points will be shown, together with a discussion of the results. 

II. KINEMATICS AND AMPLITUDES OF THE PROCESS h tH' 
A. Kinematics 

The kinematics of the process b g ^ tH~ is expressed in terms of the h quark momentum 
Pb, hehcity A^, spinor u{ph, ^h), the t quark momentum pt, helicity At, spinor u{pt, At), with: 

Pk = {Eb, 0, 0, p) pt = {Et; p' sin 0, p' cos B) , (1) 

the gluon momentum p^, helicity Ag, polarization vector and the Higgs boson momentum 

Pg = {p; 0, 0, -p) egiXg) = (0; ^, - 0) (2) 
Ph = {Eh', —p sin 9, 0, —p cos 9) (3) 
We also use the s-channel and u-channel momenta: 

q = Pg + Pb = Ph + Pt s = q' = p^ - pg = pf^ - pfj u = q"^ (4) 
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The invariant amplitude of the process b g ^ tH~ will be decomposed on a set of 8 forms 
Jkri, where rj represents the chirality R,L (sometimes denoted rj = +1,-1). The 8 scalar 
functions A''fc^(s,t,M) will be computed in the next subsection from the various Born and 
one loop diagrams. 

A = '^JknNknis,t,u) (5) 

k 

Jlr, = ig^Pr) ^2r, = {e.pt)Pr^ (6) 

Jsr, = ^Pr, Jatj = {e.pt)^gPr^ (7) 

withP^ = P^,z. = (l±7')/2. 

The 8 helicity amplitudes Fx^^XgM ci-re obtained from Dirac decompositions of the 8 in- 
variant forms. 

Averaging over initial spins and colours and summing over final spins and colours with 

col 



one gets the elementary cross section: 

da (3' 



\Fxt,Xg,Xt\'^ (9) 



spms 

where p = 2p/y^, (3' = 2p'/y^. 



B. Born and one loop amplitudes 

The Born terms result from the s-channel b exchange and the u-channel t exchange of 

Fig.m 



^Born s ^ _ ( ( [c^ (ft ^ tif-)P^ + c^(6 _ tH-)Pn]{^ + m,)^w(&) (10) 



leads to the scalar function 



and 



^Born u ^ _ ( ( |)^(^)^(^/+ ^^[^^(^ ^ ^^-)p^ ^ ^R^^^ ^ tH-)Pa]u{b) (12) 



4 



to 

N?r " = -'j4 ft^LV (13) 



with the btH~ couphngs 

c [b^tH } - — — c (6 ^ tii ) - -= — (15) 

The one loop EW terms can be classified as: 

— counter terms for b, t, H~ lines and btH~ coupling constants. We follow the on-shell 
scheme in which all of the counter terms can be computed in terms of self-energy 
diagrams. For what concerns the H~ line and the btH^ coupling, we use the procedure 
given in 



191 ] which takes into account the G , H mixing and expresses the counter 
term for tan /3 in terms ofW — H mixing self-energy. Other procedures, e.g. {20] or 21|. 



would lead to a similar divergence cancellation but differ by minor finite contributions. 

— self-energy corrections for internal b and t propagators; 

— s-channel left triangles: (Vqq), {Sqq), (xQQ) and right triangles: {qSV), {Vq'q), {qVS), 
(Sff), ifSS'); 

— u-channel up triangles: (ffV), (ffS), (SSf); and down triangles: (VSq), (q'qV), 
(SVq), if'fS), (S'Sf); 

— direct boxes: (qqq'V), (qqq'S), {qqq'xj)'i crossed boxes: {qqVS), {qqSV), (qqXiXj): 
{qqSS'); twisted boxes: {qqq'S), {qqq'V), {qqq'Xj)- 

All these contributions have been computed using the usual decomposition in terms 
of Passarino-Veltman functions and the complete amplitude has been implemented in the 
numerical code PumaMC. 

We have checked the cancellation of the UV divergences among counter terms, self- 
energies and triangles. This cancellation occurs separately inside 8 sectors, i.e. s-left L or 
R, s-right L or R, u-up L or R, u-down L or R. 
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Another useful check can be done using the high energy behaviour of the amphtudes. 



High energy rules 



22 



231] predict the logarithmic behaviour of these amplitudes at one loop 
level. They use splitting functions for external particles b, t, H and Renormalization Group 
effects on the parameters appearing in the Born terms. 



By using the logarithmic expansions of the Passarino-Veltman functions [2J] we have 
checked that the amplitudes obtained by summing the contributions of the above self-energy, 
triangle and box diagrams satisfy these rules. 

At high energy we first observe the mass suppression of N^°^"' as well as the cancella- 
tion of A^f °™ One remains with only A"!^"™ " and the 2 helicity amphtudes 

pBorn egsmtCot/3 X' e 1-cose 

-'+'+ swMw 4^'°'24 + cos^^ ^ 



-^Born 
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23| predict the corrections: 



egsmbtaxi(3 \\ 6* 1 — cos 6* . . 

At one loop logarithmic level the aforementioned rules 

= {1 + \[cihh L) + citt R)] + c^-^^H-)} (18) 

= Ff- {1 + \[c{hh R) + c{U L))] + c^-_^_{H~)} (19) 

in which cibb L), c(ti R) represent the b,t splitting functions and c'^^±^±{H~) the total 
of the H~ splitting and of the parameter renormalization of the btH~ couplings through 
6g/g — 6M/Mw + Snit/mt — 5 tan/3/ tan /3 and 6g/g — 6Mw/Mw + Smb/rrib + 5 tan/5/ tan/?. 
The result being: 



ff,-"{lH-[^|{-|^ 



a,. .l + 2c2,. 



1 2 s 
log -T 



[log' 



-t 



1 



rrir 



log^ 



w 



[log' 



-u 



}} 



(21) 



with the absence of linear logarithmic terms as noticed in 



25| 
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Taking our complete one loop computation and retaining only the logarithmic parts of 
the B,C,D Passarino-Veltman functions appearing in the various diagrams, we do recover 
the above expressions for the 2 leading amplitudes. 



C. QED radiation 



The computation of the real photon radiation contributions has been performed accord- 



ing to Ref. 
FeynArts 



26|. The matrix e 



ement has been calculated analytically with the help of 



271] and FormCalc 28|]. Infrared (IR) singularities have been regularized within 



mass regularization, i.e. giving a small mass to the photon, and the phase space integration 
has been performed using the phase space slicing method. 

Concerning the choice of the parton distribution functions (PDFs) and their factoriza- 



tion, we follow Ref. 



26j. The PD 



distribution functions CTEQ6L 



s used through this computation are the LO QCD parton 



291 ]. The factorization of the bottom PDF has been per- 



formed in the MS scheme at the scale Q = {rrit + ttih-). If the DIS factorization scheme is 
used, the differences in the numerical value of the one-loop EW effects are of the order of 
0.01% in all the considered mSUGRA benchmark points. 

The phase space slicing method introduces a fictitious separator AE in the integration 
over the photon energy. As a check of our computations we have verified that, for sufficiently 
small AE values, the final cross section is independent on the choice of AE. Despite the 
strong sensitivity to AE of the soft and of the hard cross section, c.f. the upper panel of 
Fig. [21 the dependence of the total result on AE is far below the integration uncertainties 
(lower panel of Fig. [2]). 



III. ONE LOOP RESULTS 



For the numerical evaluation of the one-loop corrections we have considered as SM inputs 
the values in TabUl The strong coupling constant has been evaluated at the renormalization 
scale Q = rrit + rrifj- and its numerical value will be given below. Since we have performed 
our computations in the on-shell scheme, we have evaluated the pole mass of the bottom 
quark starting from the MS mass at NLO in QCD, obtaining = 4.58 GeV. 
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Coupling constants 


a = 1/137.035999 as{Mz) = 0.118 


Gauge boson masses 


Mw = 80.424 GeV Mz = 91.1876 GeV 


Quark Masses 


ruu = 47 MeV = 1.55 GeV nit = 170.9 GeV 
rrid = 47 MeV m, = 0.15 GeV m^(r7l^) = 4.2 GeV 


Lepton Masses 


me = 0.51099906 MeV = 105.6583 MeV = 1.777 GeV 



Table I: Numerical values of SM inputs 



As a first step, we have analysed the distributions of the invariant mass of the final 
states da /dMinv and the total cross section for a couple of representative SUSY benchmark 
points (assuming a mSUGRA supersymmetry breaking): LS2 3^ and SPSla 31|- The 
characteristics of the benchmark points, together with the mass of the charged Higgs H~ 
and the value of as{Q)i are shown in Tab. [TTl The two benchmarks are characterized by 
largely different input parameters at GUT scale, leading to different scenarios for low energy 
spectra: the LS2 point is an optimistic "light SUSY" scenario, while the SPSla point is a 
standard and widely studied scenario for phenomenological analyses with higher masses. 
Moreover, the two points differ for the tan/3 values: LS2 features a very large tan/5 = 50, 
while in SPSla tan/3 = 10. The complete spectra at low energy have been obtained running 
the parameters through the code SUSPECT[32]. The values of tan/5 at low energy have been 



translated from those obatained in the DR scheme used by SUSPECT to the values in the 



on-shell scheme through the relation iaxi(3{0S) = tan /3{DR) — 5tan/5(OS')|finitc- The values 
we have obtained for LS2 and SPSla are 60.5 and 10.4 respectively. 



mSUGRA scenario 


mo 


mi/2 


^0 


tan f3 


sign ji 


H- 




LS2 


300 


150 


-500 


50 


+ 


229.6 


0.0965325 


SPSla 


100 


250 


-100 


10 


+ 


412.1 


0.0922963 



Table II: input parameters for the mSUGRA benchmark points and mass of the charged Higgs H 
(all values with mass dimension are in GeV) 



The resulting total cross sections and K-factors (where, as usual, K = cri-^ioop/ o'Bom) 
are shown in Tab lllll Due to the very mild dependence of our calculations on the PDF 



factorization scheme, only the results obtained in the MS scheme are shown. We have 
performed the analysis considering both the whole supersymmetric spectra (labelled "SUSY" 
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in the following discussion) and the "SUSY constrained" two-Higgs-doublet-model (2HDM) 
scenarios obtained from the original spectra considering only loops involving Higgs bosons 
and SM particles {i.e. without charginos, neutralinos and sfermions). 



mSUGRA scenario 


^Born 


SUSY 


2HDM 


o'l-ioop K-factor 


ai^ioop K-factor 


LS2 
SPSla 


5.589 
0.04207 


4.545 0.813 
0.04145 0.985 


5.867 1.050 
0.04170 0.991 



Table III: Total cross sections (in pb) at Born and loop level and K-factors 



It is possible to see that for both LS2 and SPSla the corrections in the 2HDM subset 
are very small, of the order of few percent, while in the complete SUSY case the light LS2 
spectrum features a bigger correction (~19%) than in the SPSla case (~2%). 

The differential distributions for the two benchmark points are shown in Fig. [3|, where it 
is possible to see that, as a general behaviour, the one- loop corrections decrease from the 
low invariant mass region to high energies. In the SPSla case (both SUSY and 2HDM) the 
one-loop corrections are positive near threshold, but suddenly drop and become negative 
at high energies: such compensating contributions are at the origin of the small one-loop 
correction to the total cross section. In LS2, on the other hand, SUSY and 2HDM behave in 
different ways: in the former case the one-loop corrections are always negative, the K-factor 
is ~0.97 near threshold and decreases at high energies with a behaviour analogous to the 
SPSla case, thus explaining the large negative correction to the total cross section in this 
scenario; in the latter, the one-loop corrections are positive in a wider Minv range, giving 
rise to the positive overall correction to the total cross-section. 

As a second step in our analysis, given the relevance of tan f3 for the process under inves- 
tigation, we have also analyzed the dependence of the K-factors on this parameter. Starting 
from the two previous LS2 and SPSla spectra, we have considered tan/5 as a free parameter 
and varied it at low energy within a reasonable range. Chargino and neutralino masses and 
mixing matrices depend on the value of tan/?, and they have been varied accordingly. The 
results of the K-factors as a function of tan/3 for the LS2- and SPSla-like spectra are shown 
in Fig. m It is possible to notice that the dependence of the K-factor is stronger in the LS2 
case: in the complete SUSY scenario, it ranges from ~1 (low tan/3) to ~0.89 (large tan/3), 
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while the dependence in the 2HDM scenario shows opposite behaviour. One sees that the 
NLO effects remain perturbatively under control in the whole considered range, even for 
large values of tan/3, where corrections are usually expected to become large. Similar re- 
sults have been found for the SPSla-like spectra, where however, the dependence is milder 
than in the LS2-like case. 

As a final remark, we can say that the dependence on the factorization scheme of our 
results is very mild, of the order of 0.01% in all the considered cases, and in Fig. [5] the 
differences between the two schemes are shown in more detail. 

IV. CONCLUSIONS 

In this paper we have calculated the complete EW NLO expression of the bg tH~ 
process both in a 2HDM and in the MSSM, assuming a mSUGRA symmetry breaking 
scheme, to investigate the size of the corrections to tree-level observables and their tan (3 
dependence. In our calculation we have included the full computation of QED radiation, 
which makes our analysis testable against future data. We have considered two benchmark 
points characterized by quite different values of tan/3 (10 and 50) and we have let the 
parameter vary into a reasonable range to investigate for dependences of the observables. 
We have found that the NLO corrections to the total cross sections can be sizable (negative 
and of the order of 20%) in the LS2 point, which is characterized by a light spectrum, and 
due to its cross sections (~5 pb) they might be hopefully observed at the LHC. 

The dependence on tan jS of the corrections is similar in the two benchmark points that 
we have analyzed, but more enhanced in LS2. On the other hand, the corrections exhibit a 
different behaviour in the two considered physical scenarios: in the 2HDM the corrections 
are generally mild, of the order of a few relative percent in the whole scanned range, and the 
effects raises for large tan/3; in the MSSM case, the one-loop corrections become negative 
and decreasing for large values of tan/3. 

Given the outcome of our computations, we conclude that a complete calculation of 
EW MSSM NLO effects is worth and should be taken into account for a full, reliable and 
meaningful NLO analysis of this important process, which is probably the only one that can 
provide information on the charged Higgs couplings of the model. 
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Appendix A: COUNTER TERMS AND SELF-ENERGY CORRECTIONS TO 
BORN AMPLITUDES 



In this appendix the expression of the countermterms are exphcitely hsted. They concern 
the counterterms for b, t, H~ hnes as well as the propagator self-energy corrections for h 
and t exchanges. 

s-channel counterterms 



Nit ' = - ^^il^Zic'^ib tH~) + 1(54, + #*)c"(6 - tH^ 



^bc\h ^ m-) + \Y,bZ]^c\h ^ tj)} (Ai) 



1 



+<5c^(6 ^ tE-) + - 5^5Z;iC^(6 ^ tj)} (A2) 



' = ^^^(^^^ - ^ZDc^ib - tH-) (A3) 



AT- ^ = "^^^{SZl - ^Z\y\h tH-) (A4) 
s 

where, because of the H- — G~ mixing, we denote by j = 1 and G- by j = 2. 
And from b s.e. one gets {rj = +1, —1 refering to R, L chiralities): 

+2ml{EUs) - \{6Zl + 6Z\) - ^)] (A5) 



— [E^{s) + dZ^ + Es{s) 



[s — m 



h8Z'^ + 5Zt, 



6mb 
rrib 



(A6) 
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u-channel counterterms 



l\TCt U 



2' {u-m1) ''2 



+5c^{h tH- 



3 



(A7) 



N^R 



,X 



5Z^j^)c^{h tH- 



-5c\b tH-) + \Y, SZ;,c^{b - tj)} 



N: 



ct u 
2L 



-6c' 



■-Mt + 7;5Zl)c^{b ^ tH 



tH-) + lj2^Z*,c\b^tj)} 



(A9) 



n: 



ct u 
2R 



-dc'^ib 



tH-) + l-J2^Z*,c^ib^tj)} 



tH- 



(AlO) 



N: 



ct u 
3L 



.\\mtC^{b^tH- 



'-{6Z'^ - 6Zi} 



(All) 



mtc^{b tH' 



-{5Zi-6Z'^] 



(A12) 



and from t s.e. one gets: 



tH' 



1' iu — rri^y 



n(S%(n) + 5Z%)+m?(S;(n) 



+2m^(S*5(«)--(5Z* + 5Z%) 

z rrit 



5Zl) 



(A13) 



(A14) 



NL- = ad) ^ 'V"^' Kin) + + - ^(^^^ + J - ^] (A15) 



(u — mf ) ^ 
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The c.t. appearing in the above expressions are obtained in terms of self-energies as follows: 
b and t quark 



6Zl = 6Zi ^ 6Z, = ~i:l{ml) - ml[T:i{ml) + S^(m^) + 2S',^K)] (A16) 

54 = -E^(m^) - ml[T:i{ml) + S'^K) + 2E'|K)] (A17) 

5Z'^ = 5Zl + Si(m?) - S^K) (A18) 

8^>, = -{T.l{ml) + 5Z, + m?[E2K) + S'^K) + 2S^(m?)]} (A19) 

5m, = ^i?e[SiK) + S5,K) + 2S|K)] (A20) 

5m, = ^i?e[SiK) + E^K) + 2EMm?)] (A21) 



gauge boson 



= -S'-(O) + 2-^S--(0) + ^[^!^ - ^!^] (A23) 



= ReT}^^{M^) 6Ml = ReE^^{Ml) (A24) 



Higgs boson 



We need SZ*-^ which means 5Z*^_^_ and We use the on-shell procedure of Wan 

et al 13] in which 

5Zh-h- = -^H- (P' = ml- ) (A25) 

and 

^^G-H- - (^^G+/f+ V7 - 77 lA2bj 

Couplings 



The Yukawa htH coupling leads to the c.t. 5c^ and 5c^, computed in terms of 5g, 
6mt^b, SM w (g iven above) and (5tan/3. For the latter, we have adopted the renormalization 



scheme of 



19|. 
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Sc^ _ 6g__^ 5rnt_ _ 5Mw _ (A27) 
g nit Myy tan/5 

— = — H — H ^ (A28) 

c" g TJib Mw tan p 



^ = SZ^ - -5Zf (A29) 
5' 2 



5tan/5 i?eSj^+vi/+(m 



2 ^ 

H+' 



tan/3 Miysin2/3 



(A30) 
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Figure 1: Born diagrams: s-channel bottom exchange and u-channel top exchange. 
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Figure 2: Upper panel: dependence of the 0{a) soft plus virtual and hard cross sections on the 
soft-hard separator AE. Lower panel: independence of the sum of 0{a) soft plus virtual and hard 
cross sections of the separator AE. 
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Figure 3: Differential distribution (upper panels) and partial K-factors (lower panels) in LS2 and 
SPSla. 
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K factors: LS2-like spectra 
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Figure 4: K-factor dependence on tan/3 for LS2- and SPSla-like spectra. 
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Figure 5: Factorization scheme dependence of the K-factor 
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